Ohio Biological Survey Notes 10: 1-8, 2021. © Ohio Biological Survey, Inc. 


Snake Community Composition in Unmanaged Grasslands of Different Successional Ages at the 
James H. Barrow Biological Field Station 


MEREDITH S. FITSCHEN-BROWN!, ANDREW D. RUNYON?, MATTHEW W. SorRICK?, JENNIFER M. CLARK? 


'Biological Sciences Department, Ohio University, Athens, OH 45701, U.S.A.; "Natural Resources and Environmental Sciences Department, 
University of Illinois, Urbana, IL 61801, U.S.A.; *Biology Department, Hiram College, Hiram, OH 44234, U.S.A. E-mail: clarkjm@hiram. 
edu. 


Abstract: The historical destruction and fragmentation of grassland habitat has led to negative consequences for 
vertebrate taxa that rely on these ecologically rich ecosystems. Although the decline of some taxa such as grassland 
birds has been well documented, little is known about how grassland decline impacts snake communities. This study 
was conducted to assess differences in snake abundance and species richness in unmanaged grasslands of three 
successional ages at the James H. Barrow Biological Field Station (Portage County, OH): young (4—7 years), mid- 
successional (12—15 years), and old (>50 years). Ten corrugated metal and ten plywood cover boards were used at each 
site to address community and species-specific preferences for board type. Surveys were conducted weekly from June 
through August over four summers (2016-2019) in an effort to detect seasonal peaks in abundance and the effects of 
board maturity. Additionally, snake assemblages were compared during morning and evening surveys during a single 
summer to address differences in time of day. A total of 419 eastern garter snakes, 85 northern brown snakes, 172 
eastern milk snakes, and 3 northern red-bellied snakes were captured across all years. Overall, both abundance and 
species richness were significantly higher in the youngest successional grassland, with the mid-successional and oldest 
successional site yielding similarly low diversity. Although both eastern garter snakes and eastern milk snakes were 
found in highest abundance in the youngest site, northern brown snakes were equally abundant in grasslands of each 
successional age. In the youngest successional site only, both abundance and species richness significantly increased 
across years, peaking at year three and remaining similar during year four, indicating that a threshold of diversity may 
be reached between 7 and 12 years of successional age. Because this trend was only observed in the youngest site 
and seasonal peaks were not observed for any species, it is unlikely that board maturity plays a role in detection of 
snakes. No significant differences in snake abundance or species richness were observed for cover board type or time 
of day of survey. Overall, our results suggest that both metal and plywood cover boards can be effective to survey 
snake assemblages, and evening sampling could yield similar results to morning sampling. However, these results 
are limited to one geographical location. Further, the maintenance of young successional habitat may be necessary to 
promote grassland snake diversity, but species-specific considerations should be taken into account. 
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Introduction 


Throughout Ohio’s history, its landscape has experienced dramatic shifts in the amount of cover by grassland habitat. Prior 
to European settlement, grasslands were rare, covering only 2.5% of the state. This steadily increased post-settlement with 
conversion of forest to agriculture, and has been declining over the last 50 years (Ohio Division of Wildlife, 2015). An estimated 
seven million acres of agricultural land has been lost across Ohio since 1950, with conversion to urban landscape being the 
primary cause since the 1990s (Western Reserve Land Conservancy, 2020). Subsequently, species dependent on this early 
successional habitat have been experiencing declines, with grassland birds being some of the most severely impacted taxa 
(Brennan and Kuvlesky, 2005). In response, grassland restoration practices have become common, generally with the goals 
of increasing plant diversity and improving habitat for wildlife. To achieve vegetation diversity, counteract forest succession, 
and control invasive and weedy species, a variety of management techniques have been used (e.g., controlled burns, mowing/ 
haying, grazing, and herbicides), all which can have variable direct and indirect effects on native fauna. For taxa with limited 
dispersal capabilities (such as snakes), management strategies such as mowing and fire can cause direct mortality (Durbian, 
2006), but it is suggested that the overall rate is low (Russell et al., 1999). Negative indirect effects can also occur in the 
short term, leading to unsuitable microhabitat (McLeod and Gates, 1998; Durbian, 2006) and predation risk (Durbian, 2006). 
However, some studies suggest that management may benefit snake populations or have no detectable effect. For example, 


smooth snake (Coronella austriaca) densities were similar between sites receiving brush management (mowing and tree/ 
bush cutting) versus no management (Graitson et al., 2020), and Oklahoma pastures receiving fire combined with herbicide 
treatments resulted in higher snake abundance when compared to unmanaged control sites (Jones et al., 2000). Although some 
management strategies may have minimal negative short-term impacts, the long-term effects of maintaining early successional 
habitat likely outweigh these consequences as forest encroachment leads to species replacement. 


Early successional habitat can be critical for the persistence of some snake species (Fitch, 2006a; Steen et al., 2015; Diaz and 
Blouin-Demers, 2017), but few studies have addressed how succession of unmanaged grasslands impacts snake communities. 
In northern latitudes, small colubrid snakes are found in higher abundance tn old fields than forests because old fields support 
their thermal requirements (Diaz and Blouin-Demers, 2017). A long-term study in Kansas suggested that as grasslands progress 
in successional age, there is a general trend for both increased abundance and species diversity until forest communities take 
over (Fitch, 2006b). However, these responses can be species-specific, and environmental factors may play an overarching role 
in their population dynamics (Fitch, 2006b). As agricultural lands are abandoned, these early successional areas could serve 
as important habitat for snakes. As reptiles face global declines (Gibbons et al., 2000; Reading et al., 2010), understanding 
the responses of these taxa to management and restoration strategies versus unmanaged succession is important for successful 
conservation. The goal of this study was to assess how successional age of unmanaged grasslands and cover board type impact 
snake community composition (species abundance and richness). Plywood and corrugated metal cover boards were used in 
combination to maximize sampling effort and to determine species-specific preferences, as there is controversy in the literature 
as to which type of artificial cover is most effective for sampling snakes. Further, we addressed whether community composition 
differed between morning and evening surveys (during one of the sampling years), as evening sampling is rarely reported in 
the literature. Surveys were conducted over four summers at the James H. Barrow Biological Field Station (JHBBFS; Hiram 
College, Portage County, OH, USA) in grasslands of young (4—7 years), mid-successional (12—15 years), and old (>50 years) 
successional age. We expected that the relative abundance and species richness of snakes would increase with successional 
age. Regardless of site, we predicted that snake diversity would increase from year one of the study to year four, as the boards 
became established at these sites. Further, we hypothesized that snake diversity would be higher during morning sampling 
efforts than in the evening due to morning thermoregulation behaviors. 


Methods 


Study Sites. Hiram College’s James H. Barrow Biological Field Station (Portage County, OH, USA, 41°29’90”" N, 81°10'91" 
W) covers approximately 500 acres and includes a mixture of un-timbered, old-growth beech-maple forest, transitional second- 
growth forest, successional (unmanaged) and managed grasslands, wetlands and streams. While approximately 350 acres is 
contiguous, the Eagle Creek Restoration Site property (152 acres) is separated from the rest of the field station property by 
a State route and small residential property. Two of the grasslands chosen for this study (Black Walnut Grove and Stavenger 
Field) occur on the main portion of the field station property, while the third occurs on the Eagle Creek Restoration Site property 
(Eagle Creek Grassland). 


Each grassland site was selected for surveys based on age when agriculture and subsequent years of mowing ceased or when the 
grassland was created and seeded, all remaining unmanaged since. Black Walnut Grove represents the oldest successional age 
(>50 years); Stavenger Field represents mid-successional age (12—15 years), and Eagle Creek Grassland represents the youngest 
successional age (4—7 years). The Eagle Creek Grassland was created from substrate excavated from a stream restoration 
project on the property and seeded with a native prairie mix purchased from Ohio Prairie Nursery (Portage County, Ohio, 
USA) four years prior to this study. While Solidago spp. was the dominant vegetation at all sites, the following vegetation was 
also abundant: Rosa multiflora (youngest and oldest successional sites), Rubus spp. (youngest and oldest successional sites), 
Eryngium yuccifolium (youngest successional site), and Aster vimineus (mid-successional site). 


Transect Design and Snake Surveys. A combination of 10 corrugated metal (tin roofing material) and 10 half-inch plywood 
cover boards (all 1 m7) were placed at each site during early May 2016 (a total of 20 at each site; 60 total cover boards). At 
each site, cover boards were arranged in two parallel transects (10 boards per transect) spaced two meters apart horizontally 
and five meters apart along the length of the transect. Placement of cover boards was randomly assigned to ensure 10 cover 
boards of each type were in each transect and that metal and plywood cover boards were not directly next to each other in their 
horizontal orientation. 


Snake surveys were conducted once a week during the summers of 2016 through 2019 covering a 12-week time period 
beginning the first week of June and ending in August for each year. All sites were surveyed on a single day during the hours 
of 8:30-10:30 a.m. and the order in which they were surveyed was randomly selected. They were conducted on days with 


partial to full sun and no rain during the time of survey. To survey snakes, each board was lifted, and snakes were captured by 
hand. Each individual was identified to species and measured for snout-to-vent length (cm) using a standard measuring tape. 
All snakes were released back under the cover board under which they were found. All uncaptured snakes (6%) were able to 
be identified to species. 


To compare species diversity during morning versus evening, surveys were conducted from 22 June 2016 to 30 September 
2016 following the same procedure and using the same study sites and transect design. Evening surveys took place during the 
hours of 6:30-8:30 p.m. once each week and were done on different days of the week than morning surveys to minimize snake 
disturbance. 


Statistical Analysis. A total of 172 eastern milk snakes, 85 northern brown snakes, and 419 eastern garter snakes were included 
in analyses. northern red-bellied snakes were excluded from analyses because only three were collected, and these only in the 
first year of sampling. Repeated captures were removed from analyses and were based on length and distinguishable markings 
(i.e., developmental abnormalities or injuries).Generalized linear mixed-effect models (GLMMs) using R package “Ime4” 
(version 1.1.21; Bates et al., 2015) were used to determine whether cover board type (metal versus plywood), year (2016, 2017, 
2018, 2019), or site type (youngest, mid-, and oldest successional age) influenced snake abundance and species richness. Board 
number was selected as a nested effect due to having repeated measures within a single site and a Poisson distribution was 
chosen for our inclusion of count data. A Tukey-Kramer post-hoc test was used to determine site-by-year interactions on snake 
abundance and species richness.To assess the effect of morning and evening sampling in 2016, separate GLMMs were used to 
determine differences in snake abundance and species richness, selecting site type, cover board type, and time of day as factors 
in the model. This model was constructed with board number as a random factor and a Poisson distribution was chosen for our 
inclusion of count data. A Tukey-Kramer post-hoc test was used to determine the time of day by cover board type interaction 
on snake abundance and species richness.To assess seasonal peaks from 2016 through 2019, we used a separate GLMM to 
determine if the week over the season influenced snake abundance. This model was constructed with board number as a random 
factor and a Poisson distribution was chosen for our inclusion of count data. 


Species-Specific Differences. To determine the potential association between each snake species and cover board type, separate 
chi-square tests of independence were used for each species. Data were pooled across years to determine differences in overall 
species abundance across sites. A one-way ANOVA was also used to determine species-specific associations with successional 
age, testing the difference between the abundance of each species of snake captured and site type, followed by a Tukey-Kramer 
post-hoc test to determine species by site interactions. All analyses were conducted in R version 3.4.1 (R Core Team, 2019). 


Results 


Snake Abundance and Species Richness. A total of 679 snakes (including recaptures) were recorded, representing four 
different species across all sites and years: 419 eastern garter snakes (Thamnophis sirtalis sirtalis, snout—vent-length [SVL] 
range = 13.0-93.7 cm), 85 northern brown snakes (Storeria dekayi dekayi, SVL range = 14.5—39.0 cm), 172 eastern milk 
snakes (Lampropeltis triangulum, SVL range = 22.0—98.2 cm), and 3 northern red-bellied snakes (Storeria occipitomaculata 
occipitomaculata, SVL range = 15.1—22.2 cm). Northern red-bellied snakes were captured in the mid- (two individuals) and old 
(one individual) successional sites, but not in the youngest site and only during the first year of survey. During the first year of 
capture (2016), 43 and 45 snakes were captured during morning and evening sampling, respectively. 


Both snake abundance and species richness varied across sampling years and site types. Regardless of year, the youngest 
successional site had significantly higher snake abundance (GLMM, p<0.001; Fig. la) and species richness (GLMM, both 
p<0.001; Fig. 1b) than both the mid- and oldest successional sites. Overall, species abundance was lowest in 2016 (GLMM, 
p<0.001; Fig. 1a), followed by 2017, (GLMM, p<0.001; Fig. 1a), while 2018 and 2019 did not differ in species abundance 
(GLMM, p>0.05; Fig. 1a). Similarly, species richness was significantly lower in 2016 than any other sampling year (GLMM, 
p<0.001; Fig. 1b), and 2017 had significantly lower species richness than 2018 (GLMM, p=0.023; Fig. 1b) and 2019 (GLMM, 
p<0.001; Fig. 1b). Snake abundance (GLMM, p<0.001; Fig. 1a) was significantly lower at the youngest successional site in 
2016 than all other sampling years. Species richness was significantly lower in 2016 for the youngest successional site than 
in 2017 (GLMM, p<0.001; Fig.1b). However, no significant differences in abundance (GLMM, p=0.58; Fig. 1a) or species 
richness (GLMM, p=0.71; Fig. 1b) occurred across all four sampling years for the mid-successional and oldest successional 
sites. Regardless of site type and year of survey, cover board type (metal versus plywood) did not have a significant effect 
on either snake abundance (GLMM, p=0.43; Fig. 2a) or species richness (GLMM, p=0.33; Fig. 2b). Further, neither snake 
abundance (GLMM, p=0.83; Fig. 3a) nor species richness (GLMM, p=0.99; Fig. 3b) varied based on time of day (morning 
versus evening surveys). However, in 2016, significantly more snakes and higher species richness were detected under plywood 


than metal cover boards (GLMM, p=.024, and p=0.013, respectively; Fig 3a, b). Additionally, there was no effect of week on 
snake abundance or species richness over the four years (GLMM, p>0.05). 


Species-Specific Differences. No association was detected between snake abundance and cover board type for any of the 
three main species captured (eastern garter snake, northern brown snake, and eastern milk snake; chi-square test, y7(4)=6, 
p=0.19). However, there was a significant effect of successional site age on species-specific abundance (ANOVA, F, ...,= 
73.56, p<0.001; Fig. 4a-c). Significantly more eastern garter snakes than northern brown snakes were captured at the youngest 
successional site (ANOVA, p<0.001). Additionally, eastern garter snake abundance was higher at the youngest successional 
site than either the mid- or oldest successional sites (ANOVA, both p<0.001; Fig. 4a). The youngest successional site 
also had a higher abundance of eastern garter snakes compared to eastern milk snakes (ANOVA, p<0.001; Fig. 4a and c), 
and abundance of eastern milk snakes was significantly higher at the youngest successional site compared to the oldest 
successional site (ANOVA, p=0.035; Fig. 4c). All other comparisons of species abundance across sites were not significant 
(ANOVA, all p>0.05). 
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Figure 1. Average snake abundance (a) and species richness (b) (mean + 1 SE) across all four survey years for each site type. Data are 
pooled across cover board type. Black bars = 2016, dark gray bars = 2017, charcoal bars = 2018, light gray bars = 2019. The asterisk 
indicates a significant difference by site type. 
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Figure 2. Average snake abundance (a) and species richness (b) (mean + 1 SE) across all four survey years for each board type. Data are 
pooled across site type. Black bars represent corrugated metal cover boards and gray bars represent plywood cover boards. 
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Figure 3. Average snake abundance (a) and species richness (b) (mean + 1 SE) during morning and evening surveys for summer 2016 
for each board type. Data are pooled across all three sites. Black bars represent corrugated metal cover boards and gray bars represent 
plywood cover boards. Letters above SE bars indicate post-hoc Tukey-Kramer results. 
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Figure 4. Average snake abundance (mean + 1 SE) for 
eastern garter snake (a), northern brown snake (b) 
and eastern milk snake across each survey year. Data 
are pooled across board type. Black circles = youngest 
successional site, white circles = mid-successional site, 
black upside-down triangles = oldest successional site. 


Discussion 


Across four years of survey, both snake species abundance and richness 
was highest in the youngest successional site. Both the mid- and old 
successional sites yielded similar diversity, suggesting that abundance 
and species richness may peak before 12 years of grassland successional 
age. Since snakes were not marked, we are unable to be certain that 
individuals were not recaptured during surveys. However, snakes that 
had the same body length were excluded from analysis, and our results 
should reflect comparative differences in abundance across sites. These 
results are supported by a long-term study in Kansas in which ceasing 
of grazing and cultivation led to increased abundance of the majority of 
herpetofauna initially, with declines occurring over time as the forest 
community encroached (Fitch, 2006a). For small colubrid snakes 
like common garter snakes and red-bellied snakes, open oldfields are 
preferred over forested habitats as they offer more favorable thermal 
conditions (Diaz and Blouin-Demers, 2017), further suggesting that 
as grasslands age, snake species relying on this habitat are negatively 
affected. 


In the youngest successional site, both abundance and species richness 
increased from year one of the study, peaked at year three, and then 
remained similar across years three and four, suggesting that populations 
may have been establishing in this early successional site. Since 
abundance and species richness remained similar for the mid- and old 
successional sites across the four-year study, it is unlikely that length of 
time that cover boards were present played a role in detection of snakes 
and further supports that there is likely a threshold reached in snake 
diversity that correlates to successional age. Supporting this, another 
study suggested that newly placed cover boards take only two months 
to reach peak efficiency (Grant et al., 1992). Responses to successional 
age were species-specific, both in this study and those conducted by 
Fitch (2006a, 2006b). In this study, eastern garter snakes and eastern 
milk snakes were the most common species detected. Although eastern 
garter snakes were in highest abundance at the youngest successional 
site, eastern milk snakes were found in similar abundance in the young 
and mid-successional sites. Similarly, Fitch (2006a) reported that eastern 


milk snakes increased for more than a decade and then declined with successional age. In contrast to our study, Fitch (2006b) 
reported that eastern garter snakes did not undergo drastic declines like the other snake species in their study and showed 
fluctuations in occurrence over time. The close proximity of our young successional site to aquatic habitats may be driving this 
pattern, as this species feeds largely on aquatic prey (Fitch, 2006a). However, northern brown snakes were found in similar 
abundance across all three sites, suggesting that either successional age did not drive their habitat use, or (because they were 
detected in low numbers overall) changes in population may not be able to be observed. This same pattern was observed in 
this species by Fitch (2006b), which also had very low numbers of occurrences and suggested that population changes could 
not be determined due to low sample size. This species, in particular, is cryptic and difficult to capture due to its fossorial 
behavior. northern red-bellied snakes were only detected during one year of sampling, with very few individuals captured. 
Surveys conducted at JHBBFS between 2001 and 2004 also captured all four of these species and, similar to our study, found 
eastern garter snakes in high abundance (Meshaka et al., 2008). However, in contrast, they found northern brown snakes in 
high abundance and eastern milk snakes were detected in substantially lower numbers than the former two species (Meshaka 
et al., 2008). Also, similar to our results, red-bellied snakes were present, but in low numbers compared to the other three 
species (Meshaka et al., 2008). Because individuals were not tracked in our study, it is difficult to determine whether northern 
brown snakes have declined and eastern milk snakes have increased over time. Further, Meshaka et al. (2008) detected seasonal 
peaks by age class, which was not observed in our study due to low numbers of juveniles. Additionally, Meshaka (2010) 
detected species-specific seasonal peaks in occurrence under cover boards during certain times in summer with similar species 
in Pennsylvania, and these seasonal peaks were not observed in our study. 


A variety of abiotic and biotic mechanisms are likely driving these patterns and may be strongly correlated with successional 
age. For herpetofauna in general, thermoregulatory requirements are among the main factors impacting habitat selection because 
much of their physiology and behavior is dependent on temperature (Huey, 1991; Lelievre et al., 2011), especially for snakes 
in cold climates and temperate zones (Lourdais et al., 2013). Since range in thermal requirements is species-specific, some 
thrive in open habitats while others can tolerate habitats with a wider thermal regime (Lelievre et al., 2011). Additionally, the 
availability of cover objects, both artificial and natural, likely has a strong impact on habitat selection because they offer snakes 
protection from predation as well as potentially favorable hydric and thermal microhabitats (Grant et al., 1992; Hampton, 
2007; Leliévre et al., 2010). In this study, neither species-specific nor community responses indicated preference for plywood 
or metal cover boards, with the exception of the first year of study, when comparisons were made across morning and evening 
surveys and surveys extended into September. During that first year, both abundance and species richness were highest under 
plywood cover boards but did not differ for time of day, suggesting that seasonality may play a role in board preference type, as 
no board preference was detected during morning-only surveys from June through August. This contradicts other studies that 
show preferences by reptiles for tin cover objects over plywood (Grant et al., 1992; Hampton, 2007). Other studies suggest that 
responses can be species-specific—e.g., red-bellied snakes can show a preference for tin (Halliday and Blouin-Demers, 2015; 
Diaz and Blouin-Demers, 2017), while garter snakes have conflicting preferences in the literature (Englestoft and Ovaska, 
2000; Diaz and Blouin-Demers, 2017)—but no species-specific preferences were observed in this study. 


Although the majority of studies survey snakes in the morning, our results suggest that evening surveys could produce a similar 
outcome for population and community studies. Grant et al. (1992) found time-of-day surveys to be site-specific, also making 
concessions that a variety of factors could influence whether or not time of day had any impact on the community assemblage. 


Although temperature and cover objects are often accepted as the driving mechanisms for habitat selection, factors including 
but not limited to environmental moisture (Fitch, 2006b), food availability (Huey, 1991; Fitch, 2006b), predation risk (Webb 
and Whiting, 2005), vegetation cover (Pringle et al., 2003), and interspecific competition (Luiselli, 2006) can impact habitat use 
by snakes at both the individual and species level and are likely not mutually exclusive. For example, Fitch (2006b) suggested 
that low moisture in drought years negatively impacted the availability of earthworms (a major food source for young-of-year 
red-sided garter snakes), leading to low survivorship and possibly delayed maturity (Fitch, 2006b). Further, snake declines were 
detected in species specializing on small rodents (e.g., voles), which correlated with vole declines associated with increased 
grassland successional age (Fitch, 2006a). In our study, up to nine individuals were captured under a single board, and it was not 
uncommon to find eastern garter snakes, northern brown snakes, and eastern milk snakes sharing boards. Of particular interest, 
snakes appeared to avoid boards 1n which ants developed colonies, suggesting that a variety of factors impact habitat use by 
snakes and microhabitat-level constraints are important to consider. 


Many studies have linked landscape-level factors such as surrounding land use type (Cagle, 2008), distance from habitat edge 
(Cagle, 2008), and habitat fragmentation (Breininger et al., 2012) to snake occurrence, but responses are often species-specific. 
Further, roads also fragment landscapes and associated snake mortalities are also a known risk for populations (Breininger et 
al., 2012; Jochimsen et al., 2014). In this study, landscape features could impact source populations, mobility, and dispersal 
as a mix of forest, agriculture, and residential areas all occur within a kilometer of each survey site. Patch size can also be an 


important driver of snake community composition, with larger patches having highest abundance and species richness (Kjoss 
and Litvaitis, 2001). However, this was likely not a factor in this study, as the mid-successional site was the largest grassland 
patch and the youngest and oldest sites were both relatively small and of similar size. 


Although there is limited published literature describing the effects of successional age of unmanaged grasslands on snake 
assemblages, it is generally understood that open and early successional habitats are important for snake populations. Habitat 
loss 1s a major factor in the decline of snake populations (Gibbons et al., 2000), and with grasslands being the most endangered 
ecosystem in North America (Samson et al., 2004), more attention should be focused on taxa that are reliant on these habitats. 
Population declines are particularly evident in areas with dwindling habitat such as Midwest tallgrass prairies (Cagle, 2008). 
Management of grasslands to prevent encroachment of woody vegetation is likely necessary to provide this critical open habitat 
for some species (Lagory et al., 2009). However, several studies have suggested that these practices should be timed with 
periods when snakes are inactive to minimize mortality, and that management strategies create high habitat heterogeneity to 
benefit multiple species. Because our results are limited to only one study area, generalized conclusions cannot be made and 
more studies should be conducted to look at general patterns of successional age on grassland snake assemblages, especially 
with comparison made in unmanaged versus managed sites using long-term surveys. 
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A New Record of Burrowing Mayfly, Anthopotamus neglectus neglectus (Traver, 1935) 
(Ephemeroptera: Potamanthidae), from Ohio, USA 
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Abstract: A new state record for a mayfly (Ephemeroptera) was collected on the Olentangy River, Delaware County, 
Ohio, USA. Anthopotamus neglectus Traver (1935) were collected as nymphs and subsequently reared to adults. 


Keywords: Olentangy River, Delaware County, Ohio 


Introduction 


The neglected hackle-gilled burrowing mayfly, or the golden (or yellow) drake to fly fishers, Anthopotamus neglectus was 
first described by Traver (1935) as Potomanthus neglectus. Bae and McCafferty (1991) reorganized the family Potomanthidae 
and placed the taxon in a new genus, Anthopotamus McCafferty and Bae (1990). They further divided the species into two 
subspecies, A. neglectus neglectus and A. neglectus disjunctus. The geographic range of the former species was originally given 
as a small circle centered in New York. The latter species was centered in the south-central United States. More recently, A. 
neglectus neglectus has been reported in eastern North America including Ontario, Alabama, Arkansas, Maryland, Missouri, 
Mississippi, New York, Oklahoma, Tennessee, Virginia, and West Virginia (Randolph, 2002). The online database NatureServe 
Explorer (2019) lists the range of A. neglectus neglectus as previously stated, with the addition of Georgia and Pennsylvania 
(but it includes the caveat “Distribution data for U.S. states and Canadian provinces is known to be incomplete or has not 
been reviewed for this taxon”). It is notable that the Upper Midwest states, including Ohio, are not included in the published 
geographic distribution. The published geographic distribution of Arthopotamus neglectus disjunctus includes Oklahoma, 
Arkansas, Mississippi, and Tennessee (Randolph, 2002). 


Methods 


Anthopotamus neglectus nymphs were collected on 27-VII-2019 and 05-VHI-2019 from the Olentangy River, Delaware Co., 
Ohio, USA at Highbanks Metro Park (40.149145, -83.042546). At the sampling site, the stream bed is sand and cobble. The 
nymphs were most prevalent in the rapids. Water temperature was 20.6°C (69°F) at the time of collections. 


Sampling Methods. A standard one m? kick net was used to collect nymphs in riffles from the above sites. Mature larvae were 
placed in a bucket of stream water oxygenated with a battery-operated bubbler. Mature larvae were reared in an aquarium until 
emersion. Subimagoes emerged over a two-week period beginning 08-VIII-2019; imagoes (male and female) emerged a day 
later than the subimagoes. The imagoes lived only 2 days. Specimens were deposited at The Ohio State University Museum of 
Biological Diversity and Mayfly Central, Purdue University. 


Results and Discussion 


Our collections of A. neglectus represent the first record of this mayfly species from Ohio. Nymphs were identified to genus 
with the keys and description of Merritt, Cummins, and Berg (2008) and to species with the description and keys of Bae 
and McCafferty (1991). The identification was confirmed by Luke Jacobus (Indiana University—Purdue University [[UPUC], 
Columbus, IN). 


Nymphs. Anthopotamus neglectus nymphs are recognized as burrowing mayflies by the forward-thrusting upturned tusks 
(mandibles). The abdominal gills are oriented laterally, which distinguishes them as the family Potamanthidae (see Fig. 1). The 


particular shape of the tusks (1.e., narrowing abruptly and rather straight; Fig. 2a); the elongate narrow spine on the apical, inner 
margin of the foretibiae (Fig. 2b); and the extensively marked tergal segments (Fig. 1), identify it as A. neglectus. 


For the nymphs, the division of the species into subspecies is partly on the basis of regional geographic location. A. neglectus 
neglectus is found in the northeastern U.S.; A. neglectus disjunctus is found further south (Tennessee, Arkansas, Georgia, 
Mississipp1). 


Figure 1. Anthopotamus neglectus nymph, dorsal view. 


Figure 2. A) Tusk or mandible of A. neglectus, dorsal view; B) foretibia showing apical spine on inner margin 
(arrow). 


Adults. A more exact description of the subspecies is given for the adults (Fig 3). In general terms, the body of the adult is pale 
yellow with three tails; the middle tail (terminal filament) is shorter than the outer tails (cerci). The wings are transparent, with 
clearly marked wing veins. The subimago (dun) and imago (spinner) are similar in the coloration and clarity of the wings. The 
MP1 vein of the hindwing of A. neglectus neglectus is symmetrically forked with the MP2 vein (Fig. 4). The crossveins of both 
the fore- and hindwing are strongly infuscated (darkened with a brown tinge). The other subspecies, A. neglectus disjunctus, 
has the MP1 of the hindwing more commonly connected with the CuA vein, forming an MP-cell (Bae and McCafferty, 1991). 
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Figure 3. Anthopotamus neglectus female. 


Figure 4. Hindwing of male Anthopotamus neglectus neglectus. 


The collection of A. neglectus neglectus in Ohio is a new range extension for the taxon. A possible reason that it has heretofore 
not been reported in Ohio is because the Ohio Environmental Protection Agency (OEPA) and other groups involved in 
monitoring water quality do not carry the taxonomy beyond the level of genus (OEPA MacroinvertebrateTaxa List, 2019). The 
finding of A. neg/ectus in the Olentangy River, a habitat that is not unusual for Ohio but which has a high abundance of aquatic 
insects, including stoneflies (DeWalt, et al., 2012), might mean that the taxon has already been collected but not recognized by 
OFPA or other groups involved in water quality monitoring. Fly fishers report seeing the “yellow drake” in midsummer (Brian 
Flechsig, pers. comm. ). 


The family Potamanthidae consists of one genus in Ohio, Anthopotamus. There are four species in this genus in North America 
(Randolph, 2002). Anthopotamus distinctus and A. myops were previously reported in Ohio by Randolph and McCafferty 
(1998). Anthopotamus verticis has recently been recovered from a collection dated 1936 from Lake County (Bolton et al., 
2019). The finding of A. neglectus completes the list of taxa for Ohio. 
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Abstract: The pathogenic fungus Batrachochytrium dendrobatidis (Bd) is a leading cause of mortality among 
amphibians worldwide. To provide a large-scale assessment of Bd in Ohio (USA), we analyzed aggregated data on 
Bd infection prevalence, microhabitat, and seasonality from across the state. Skin swabs (n = 2,200) were collected 
from 26 species of amphibians (and one hybrid form) across 22 counties in Ohio between 2006 and 2020. Of these, 20 
species tested positive for Bd and at least one positive sample was found in 18 of 22 counties. Overall Bd prevalence 
among the tested samples was 17.7% (390/2,200). Frogs were infected at a significantly higher frequency (23.3%; 
283/1,212) than salamanders (10.7%; 106/988). Amphibians sampled from aquatic microhabitats were infected at 
a significantly higher frequency (frogs: 24.9%; 237/953; salamanders: 15.0%; 105/706) than those sampled from 
terrestrial microhabitats (frogs: 17.4%; 45/259; salamanders: 0.4%; 1/282). Seasonal infection frequency differed 
among species, with some species showing seasonal infection patterns and others not. All skin swabs tested for another 
emerging fungal pathogen of amphibians (Batrachochytrium salamandrivorans, n = 186) tested negative. While 
these data suggest that Bd is widespread in Ohio and has strong taxonomic and microhabitat associations, they also 
underscore how much work remains to be done, with many species still undersampled (or not sampled at all) for this 
emerging pathogen. 


Keywords: Batrachochytrium dendrobatidis, Batrachochytrium salamandrivorans, disease surveillance, frog, salamander, 
wildlife diseases. 


Introduction 


Amphibians are one of the most threatened groups of organisms on the planet, with nearly one-third of species thought to 
be at risk of extinction (Stuart et al., 2004). While the causes of amphibian declines are complex and often include multiple 
stressors (Blaustein and Kiesecker, 2002; Collins and Crump, 2009; Rohr and Palmer, 2013), an important driver of these 
losses is the emerging amphibian fungal pathogen Batrachochytrium dendrobatidis (hereafter Bd) (Berger et al. 1998; Scheele 
et al. 2017; Skerratt et al. 2007; Vredenburg et al. 2010). Overwhelming evidence indicates that this pathogen has played a 
major role in numerous extinctions and local extirpations, resulting in biodiversity loss (Collins and Crump, 2009; Skerratt et 
al., 2007; Vredenburg et al., 2010; Cheng et al., 2011). More recently, a newly discovered fungal pathogen, Batrachochytrium 
salamandrivorans (hereafter Bsal), has also been implicated in declines in some salamanders in Europe but has yet to be 
reported from the New World (Gray et al., 2015; Yap et al., 2015; Waddle et al., 2020). 


Patterns of Bd infection, including taxonomic associations (O’ Hanlon et al., 2018), environmental influences (Karvemo et al., 
2018), temporal (Watters et al., 2019) and geographic trends (Olson et al., 2013), and genetic divergence (Byrne et al., 2019), 
continue to be examined in detail by researchers. Global mapping efforts have been underway for some time (Olson et al., 
2013; James et al., 2015), but patterns of Bd infection at regional spatial scales are still incomplete. Many studies that examine 
infection status of free-living individuals are based on information from one or a small number of sites. These studies, while 
important, are often hampered by small sample sizes. To effectively monitor wildlife pathogens and to rigorously investigate 
factors influencing infection patterns, large-scale studies with robust sample sizes are needed. 


While there is often great interest in areas where wildlife disease outbreaks are occurring (epizootics; Lips et al., 2006), there is 
also a need for studies in areas where pathogens have long been established (enzootic stage; James et al., 2015; Scheele et al., 
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2017). For example, Bd has been present in the midwestern United States since at least the late 19" century (Talley et al., 2015). 
Recent evidence suggests that partial population recovery from Bd in some taxa is possible even over relatively short time 
scales (Voyles et al., 2018). Therefore, intensive studies in areas where Bd has been present for long time periods are critically 
important to examine the potential mechanisms of recovery from declines and to assess continuing changes. Here, we present 
new data and aggregate scattered published and unpublished data on the distribution, prevalence, taxonomic associations, and 
seasonal and microhabitat infection patterns of Bd infection for amphibians in the state of Ohio (USA). These aggregated data 
allow us to overcome some of the limitations of small sample sizes of individual species to search for overall patterns. We also 
assess the occurrence of Bsal in Ohio. 


Materials and Methods 


Sample Collection. To collect samples to assess pathogen infection status, we hand-captured adult amphibians in the field, 
using new sterile gloves for each individual. In aquatic environments, we used dip nets. Fewer than ten larvae were captured; 
we pool these data with those from adults because they represent a very small fraction of the total. Captured individuals were 
swabbed either for 30 seconds or 25 times each, using a sterile swab. Skin swabs were collected between 2006 and 2020 in 22 
counties in Ohio (USA) and stored in a freezer until they were processed in the laboratory (see below). Species identification 
was based on Pfingsten et al. (2013), but a small number of small P/ethodon were determined to be hybrids between Plethodon 
cinereus and P. electromorphus based on coloration and morphological criteria in Lehtinen et al. (2016). We worked at 110 
study sites across the state and, in general, all sites were sampled for amphibians only once. However, at Wooster Memorial 
Park (Wayne Co.), we repeatedly sampled amphibian populations from 2016 to 2020. This time series will be the subject of a 
separate manuscript, but the aggregate data are included here. For swabs from this site, each was assumed to be independent 
across years. In addition to these new data, our analysis also includes previously published data from Steiner and Lehtinen 
(2008), Krynak et al. (2012), Korfel and Hetherington (2014), and Sonn et al. (2019). Because each of the researchers was 
working independently of others at the time that skin swabs were collected in the field, sampling was haphazard. 


Infection quantification using PCR. DNA was extracted from the swabs using Qiagen DNEasy Blood & Tissue kits. Isolated 
DNA was tested for Bd using either quantitative real-time PCR (n = 1,077), following a modified protocol based on Boyle et al. 
(2004) or conventional PCR (n= 1,124) following Annis et al. (2004). Because of differing protocols among researchers, some 
samples were run in triplicate while others were run in singlicate. For samples run in triplicate, they were deemed Bd positive 
if two or more PCR reactions for that sample were positive. All skin swab samples were run to assess Bd infection status; a 
smaller subset (all from Wooster Memorial Park in Wayne county) were also run for Bsal. Unfortunately, because different labs 
used different standards in their qPCR reactions, infection load numbers for positive samples cannot be meaningfully compared 
among researchers. For this reason, we present only infection frequency data. 


Statistical Analysis. Our aggregate analyses pool infection data from many different species, allowing much larger sample 
sizes than would otherwise be available. However, due to uneven sampling among species, these pooled results should 
be interpreted cautiously. Multiple logistic regression was used to determine the influence of taxonomic order (frogs or 
salamanders), microhabitat type (terrestrial or aquatic), and season (spring, March—May; summer, June-August; or autumn, 
September—November) on Bd infection likelihood. Reference categories for these predictors were assigned: e.g., frogs, aquatic 
microhabitat, and spring, respectively. All predictor variables were coded as categorical and entered into the logistic regression 
model simultaneously. We also used x? goodness-of-fit tests to test for differences in infection frequency across different seasons 
and microhabitats and between order, family, and species when sample sizes allowed. In addition to aggregated analyses, we 
also used x? goodness-of-fit tests to assess seasonality and microhabitat influences for individual species that had at least n = 
30 for each season or microhabitat. Cramér’s V was used to assess the strength of association between variables. All statistical 
analyses were performed in SPSS version 27.0. 


Results 
A total of 2,200 skin swab samples from 26 species (and one hybrid form) in eight amphibian families were collected and 
processed from 22 Ohio counties. A total of 390 (17.7%) of these swabs tested positive for Bd, including swabs from 18 of the 
22 counties sampled (Table 1). Twenty of 26 species tested positive for Bd in at least one individual (Tables 2, 3), however, all 


swabs tested for Bsal (n = 186) were negative. 


As a group, frogs had a significantly higher frequency of infection (23.3%) than salamanders (10.7%; y7 = 60.893, df= 1, n= 
2,200, p < 0.001, Cramér’s V = -0.165; Table 2). Within frogs, infection frequency was not equally distributed among families 


14 


(x? = 41.494, df = 2, n = 1,234, p < 0.001, Cramér’s 
V = 0.183), with Ranidae having the highest frequency 
of infection (31.0%) and Bufonidae the lowest (11.2%, 
Table 2). Similarly, within salamanders, infection 
frequency was also not equally distributed among 
families (y? = 25.681, df=4, n=988, p<0.001, Cramér’s 
V = 0.161), with Cryptobranchidae having the highest 
frequency of infection (29.2%) and Plethodontidae the 
lowest (8.6%; Table 2). 


At the species level, Bd infection frequency was also not 
equally distributed among frog species (y? = 91.411, df 
= 10, n = 1,234, p < 0.001, Crameér’s V = 0.272; Table 
3). Lithobates catesbeianus (60.3%) and L. clamitans 
(30.3%) had particularly high infection frequencies. In 
one frog species (Hyla versicolor), no infection was 
detected. Bd infection frequency was also not equally 
distributed among salamander species (y? = 60.058, df 
= 15,n=991, p< 0.001, Cramér’s V = 0.246; Table 3). 
Ambystoma jeffersonianum (34.4%) and Cryptobranchus 
alleganiensis (29.2%) had particularly high infection 
frequency. In four salamander species (two Ambystoma 
and two Plethodon), no infection was detected (Table 3). 


Microhabitat use significantly influenced infection status. 
Amphibians in aquatic microhabitats were infected 
with Bd at significantly higher frequency than those in 
terrestrial microhabitats (7° = 64.460, df= 1, n=2,200, p 
< 0.001, Cramér’s V = 0.170). This also held true when 
frogs and salamanders were analyzed separately (frogs: 
y° = 14.307, df = 1, n = 1,234, p < 0.001, Cramér’s V 
= 0.108; salamanders: y? = 45.124, df = 1, n = 991, p 
< 0.001, Cramér’s V = 0.213; Table 4). Notably, only 
two species of frogs sampled in a terrestrial microhabitat 
were Bd positive (P. crucifer: 1/2; A. americanus: 
12/115). Similarly, of the 282 salamander samples taken 
from terrestrial microhabitats, only one was positive (E. 
bislineata: 1/60). Blanchard’s cricket frog was the only 
frog with a large enough sample size in both aquatic and 
terrestrial habitats to permit a within-species test. This 
species did not have significantly different infection 
frequencies between habitats (11/50 [22.0%] in terrestrial 
habitats and 74/393 [18.8%] in aquatic habitats; y? = 0.288, 
df = 1, n = 443, p = 0.592). Of salamanders with large 
enough sample sizes for a within-species test, infection 
was significantly more frequent in aquatic habitats in both 
species: the infection frequency in E. bislineata was 1/60 
(1.7%) in terrestrial habitats and 47/315 (14.9%) in aquatic 
habitats (77 = 6.678, df= 1, n =423, p= 0.01, Cramér’s V 
= 0.126), and in P. cinereus, infection occurred in 0/115 
(0.0%) of individuals in terrestrial habitats and 10/53 
(18.9%) in aquatic habitats (y? = 19.341, df= 1, n=178, p 
< 0.001, Cramér’s V = 0.330). 


Bd infection frequency in all amphibians was not equally 
distributed between seasons (vy? = 9.181, df = 2, n = 
1,662, p = 0.010, Cramér’s V = 0.074). For all amphibian 


Infection 
County 
% Bd Positive 
Auglaize 
Belmont 
Clinton 


Columbiana 


Cuyahoga 


Delaware 
Franklin 
Geauga 
Greene 

Hardin 

Harrison 

Jefferson 

Lake 
Marion 
Medina 
Monroe 
Preble 
Stark 
Summit 
Union 
Warren 
Wayne 
Total 


Table 1. Bd infection frequency across 22 Ohio counties in frogs and 
salamanders (n = 2,200, all species pooled). 


Infection 


oe %B 
Positive ees 
Bufonidae 
Hylidae 
Ranidae 


Total 


Ambystomatidae 


Cryptobranchidae 


Plethodontidae 
Salamanders 
Proteidae 


Salamandridae 


Total 
Table 2. Summary of frog and salamander Bd infection frequency in 


Ohio by family (n = 2,200). Infection was significantly higher in frogs 
than salamanders and not equally distributed across families. 
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Table 3. Summary of Bd infection frequency in Ohio by species (n = 2,200). Plethodon hybrids were P. 
cinereus x P. electromorphus. 


Infection 
Species 
Positive % Bd Positive 

Acris blanchardi 85 1933 

Anaxyrus americanus a2 

Hyla chrysoscelis 26.7 

Hyla versicolor 0.0 

Lithobates catesbeianus 60.3 

Lithobates clamitans 30.3 

Lithobates palustris 273 

Lithobates pipiens TA, 
Lithobates sylvaticus 
Pseudacris crucifer 
Pseudacris triseriata 
Total 
Ambystoma jeffersonianum 
Ambystoma maculatum 
Ambystoma opacum 
Ambystoma texanum 
Cryptobranchus alleganiensis 
Desmognathus fuscus 
Desmognathus ochrophaeus 


Eurycea bislineata 


Salamanders Eurycea longicauda 


Necturus maculosus 
Notophthalmus viridescens 
Plethodon cinereus 
Plethodon electromorphus 
Plethodon glutinosus 
Plethodon spp. hybrid 


Pseudotriton ruber 


Total 


species, there was a slightly higher infection frequency during the spring than in the summer or autumn, but patterns differed 
between frogs and salamanders (Table 5). Infection frequency was distributed relatively equally among seasons in frogs, with 
the highest frequency in the spring (27.0%), followed by lower (but still high) infection frequencies of 21.2% in the summer 
and 22.3% in the autumn. Only two species of frogs had large enough sample sizes in all three seasons for statistical analysis. Of 
these, Acris blanchardi did not have significantly different Bd infection frequency among seasons (v7 = 5.566, df = 2, n = 443, 
p =0.062, Cramér’s V = 0.112: Table 6), but LZ. clamitans did (y? = 8.569, df = 2, n= 409, p=0.014, Cramér’s V = 0.131; Table 
6). Overall salamander infection frequency was lower; it was also highest in spring (14.0%) and intermediate in the summer 
(11.9%), but markedly lower in the autumn (7.5%, Table 5). Only one salamander species had large enough sample sizes in 
all three seasons for statistical analysis. Eurycea bislineata did not have significantly different Bd infection frequency among 
seasons (7 = 2.338, df =2, n= 423, p =0.311, Cramér’s V = 0.074; Table 7). 
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Logistic regression of the infection data for all species pooled 
identified several important predictors of infection status (Omnibus 
y-= 138.377, df = 4, n = 2,200, p < 0.001; Table 8), explaining 
between 6.0% (Cox and Snell 1?) and 10.0% (Nagelkerke 1’) of 
the variance in infection prevalence. The strongest predictor of 
infection status was microhabitat. Individuals sampled from aquatic 
microhabitats were more than four times as likely to be infected with 
Bd compared to individuals sampled from terrestrial microhabitats 
(Table 8). Taxonomic differences also contributed significantly to 
the model, with frogs more than twice as likely to be infected with 
Bd compared to salamanders. The season in which an individual 
was sampled was also a significant predictor of Bd infection status 
(Table 8), with summer-sampled individuals 1.5 times less likely to 
be infected than spring-sampled individuals. Spring and fall were 


Infection 
Microhabitat % Bd 


Positive 

Aquatic 
Frogs Terrestrial 
Total 
Aquatic 
Salamanders Terrestrial 


Total 
Table 4. Aggregate infection frequency of frogs and 


salamanders in aquatic and terrestrial microhabitats (n = 
2,200). 


not different from one another (Table 8). 


Discussion Infection 


Attempts to synthesize what we know about Bd and its effects on its Bcd 
hosts have been crucial to identify information gaps, catalyze further 
research, and enable conservation efforts (James et al., 2015; Lips, 
2016). Our relatively large sample sizes from aggregated data allow 
for firmer conclusions than is often the case about the distribution Autumn 
and frequency of Bd infections, as well as identifying key factors Total 
influencing infection patterns. Our results indicate that Bd occurs 
in most sampled species (20 of 26) and is widely distributed across 
Ohio; infected amphibians were found in 18 of 22 sampled counties, 
but counties without Bd detection often had small sample sizes and 
further sampling would likely detect it (Table 1). Bd 1s also known Total 
to occur in other Ohio counties not directly assessed in this paper 
(e.g., Butler, Clermont, and Hamilton counties; Rumschlag and 
Boone, 2020) and probably occurs statewide. 


Spring 


Summer 


Spring 


Summer 
Salamanders 
Autumn 


Table 5. Aggregated seasonal Bd infection data for frogs and 
salamanders from Ohio (n = 2,200). 


Spring Summer Autumn 


Total Positive Total Positive Total Positive 

Acris blanchardi 31/117 26.5 34/199 17.1 20/125 
Anaxyrus americanus 8/78 10.3 6/40 15:0 0/7 
Hyla chrysoscelis 4/12 33.3 0/0 0.0 0/0 
Hyla versicolor 0/1 0.0 0/23 0.0 0/1 

19/26 73.1 6/19 31.6 10/13 

Lithobates clamitans* 27/67 40.3 57/174 32.8 40/168 
Lithobates palustris 3/5 60 0/6 0.0 0/0 
Lithobates pipiens 1/2 0/11 0.0 0/1 
14/65 eS 6.7 1/5 


Species 


Lithobates catesbeianus 


Lithobates sylvaticus 
Pseudacris crucifer 0/25 ; 0/1 0.0 1/2 
Pseudacris triseriata 0/0 : 0/0 0.0 0/1 


Total 107/398 104/491 72/323 


Table 6. Bd infection frequency in frogs by species among seasons. % Bd Positive indicates the percentage of swabs from that species 
during that season that tested positive for Bd. Bolded species had large enough sample sizes in all seasons (n > 30) to permit statistical 
analysis. Asterisks indicate species with statistically significant differences in Bd infection among seasons. 
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Spring Summer Autumn 


Total Positive Total Positive Total Positive 

Ambystoma jeffersonianum 11/31 35,5 
Ambystoma maculatum 9/53 17.0 7 3 0.0 0 
Ambystoma opacum 0/0 0.0 0/5 0.0 0/0 
Ambystoma texanum 0/30 0.0 0/0 0.0 0/0 
Cryptobranchus alleganiensis 0/0 0.0 4/21 19.0 3/3 
Desmognathus fuscus 2/12 16.7 3/34 8.8 0/10 
Desmognathus ochrophaeus 0/1 0.0 3/21 14.3 0/0 

Eurycea bislineata 4/65 6.2 20/150 13.3 24/208 
Eurycea longicauda 0/1 0.0 0/11 0.0 1/1 
Necturus maculosus 0/0 0.0 2/8 25.0 0/0 
Notophthalmus viridescens 8/20 40.0 0/9 0.0 0/2 
Plethodon cinereus 8/67 11.9 2/16 125 0/95 
Plethodon electromorphus 0/17 0.0 0/1 0.0 0/16 
Plethodon glutinosus 0/6 0.0 0/17 0.0 0/12 
Plethodon spp. hybrid 0/4 0.0 0/0 0.0 0/13 
Pseudotriton ruber 0/0 0.0 3/14 j 0/2 

Total 42/307 ; 37/311 : 28/373 


Table 7. Bd infection frequency in salamanders by species across seasons. % Bd Positive indicates the percentage of swabs from that 
species during that season that tested positive for Bd. Bolded species had large enough sample sizes (n > 30) in all seasons to permit 
statistical analysis. No species varied significantly in Bd infection among seasons. 


While Bd infection prevalence was generally low to moderate, infection was most frequent in true frogs (Ranidae) and in 
salamanders with aquatic life stages. These more aquatic amphibians are presumably exposed to Bd zoospores more frequently 
than terrestrial amphibians, resulting in overall levels of infection over four times higher in amphibians sampled from aquatic 
microhabitats than from terrestrial ones. Bd zoospores must remain moist to cause infection and disease (Berger et al., 2005; 
Johnson and Speare, 2005; Van Rooj et al., 2012), and in other regions, similar associations of Bd with aquatic microhabitats 
have been identified (Lips et al., 2006; Kriger and Hero, 2007; Greenberg et al., 2017). Infection frequency of salamanders 
was generally low across Ohio when compared to frogs (9.3% of 805 sampled individuals). However, salamander infection 
frequency in Ohio was much higher than in Tennessee and Virginia, which showed minimal infection (0.7%) in Plethodon 
salamanders (Muletz et al., 2011). Our results were more similar to those reported by Jongsma et al. (2019) from Canada 
(maximum of 12.9% in Plethodon cinereus). 


Overall seasonal infection frequency patterns in frogs eae PENS pegs 


and salamanders showed similar peaks of infection in Microhabitat 1.498 | . 47.649 ea < 0.001 4.471 


the spring, though fall infection frequency dropped 


more strongly in salamanders than in frogs (Table 5). Season 
However, our aggregated analyses may oversimplify ae 034 
the situation; each species may or may not conform 


to these general patterns, and species with large summer = .454 
sample sizes may bias the results toward patterns Order 0.824 | .124 | 44.085 aria <0.001 | 2.294 
they exhibit. For example, when looking at seasonal 

Bd infection patterns in the green frog (L. clamitans), Constant -3.529 | .241 | 213.527 al <0.001 | .0 

we found significant differences among seasons 


with lower frequency in the summer and autumn Table 8. Multiple logistic regression results of variables predicting likelihood 
compared to spring (see also Korfel and Hetherington of Bd infection, all species pooled (Omnibus Chi-square = 138.377, dies 4,n 


2014). Other studies have reported similar decreased = 2,200, p < 0.001). Significant predictors are in bold. Salamanders, spring, 
Bd infection frequency during the warmer or drier and aquatic microhabitat were used as reference categories. Odds indicates 
months of the year (Berger et al., 2004; Kriger and the odds of infection in the reference category compared to the category 


bei d (n = 2,200). 
Hero, 2007; Padgett-Flohr, 2008; Lannoo et al., ee 
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2011; Chatfield et al., 2012; Watters et al., 2019), possibly due to clearance of infections at higher temperatures that are less 
optimal for Bd growth (Longcore et al., 1999; Weinstein, 2009). Similarly, Blanchard’s cricket frog (Acris blanchardi) tended 
to have higher infection in spring than in other seasons (as documented in Sonn et al., 2019; Sonn et al., 2020) even though 
this pattern did not quite reach statistical significance in our analysis. However, the northern two-lined salamander (Eurycea 
bislineata) showed no statistically significant difference in infection frequency among seasons in Ohio, with large samples 
from all seasons (Tables 6 and 7). Thus, while some species may clear Bd infections under warmer conditions, others seem to 
maintain their infections throughout the active season. 


North America is a hotspot for amphibian biodiversity, especially for salamanders; it accounts for almost half of the world’s 
recognized species (Yap et al., 2015). Studies that monitor amphibian populations and their pathogens play an essential role 
in providing baseline infection data that document infection patterns. Aggregating data among researchers provides stronger 
baseline information and allows for a greater understanding of the interactions between amphibians and their pathogens. 
Despite the large amount of aggregated data presented here, however, only a few species have robust samples (Tables 6 and 7). 
Many species in Ohio are still sparsely sampled or have yet to be assessed for Bd infection at all. This underscores the need for 
ongoing disease surveillance to help mitigate current and future threats to wildlife populations from emerging diseases. 
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